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Using the density matrix renormalization group, we calculate the energy of a domain wall in the 
2D t-J model as a function of the linear hole density pi, as well as the interaction energy between 
walls, for J/t = 0.35. Based on these results, we conclude that the ground state always has domain 
walls for dopings < x < 0.3. For x < 0.125, the system has (1,0) domain walls with pg ~ 0.5, while 
for 0.125 < x ~ 0.17, the system has a possibly phase-separated mixture of walls with pi ~ 0.5 and 
pi = 1. For x £. 0.17, there are only walls with pi = 1. For pi = 1, diagonal (1,1) domain walls have 
very nearly the same energy as (1,0) domain walls. 



PACS Numbers: 74.20. Mn, 71.10.Fd, 71.10.Pm 

In the last few years experimental evidence for stripe 
formation in the cuprates has been mounting [EJ. A 
proper theoretical description of domain walls and 
striped phases in a doped two-dimensional antiferromag- 
net has been extremely difficult to develop, however. Al- 
though simple mean field theories for the t-J or Hubbard 
models yield domain walls, it is clear that real domain 
walls have a much more subtly correlated ground state. 
Partially filled domain walls have been particularly hard 
to describe theoretically. 

Recently, we reported numerical results [^) showing 
a striped phase in a 16 x 8 t-J system at a filling of 
x = 0.125, which were in agreement with neutron scatter- 
ing results for Lai.e-zNdo^Sr^CuC^, a system in which a 
suppression of superconductivity occurs near x = 0.125. 
Here, domain walls with a linear filling of 1/2 hole per 
unit length separating 7r-phase-shifted antiferromagnetic 
regions were spaced four lattice spacings apart. The ques- 
tion of what happens at other fillings was not addressed, 
and the possibility of other types of domain walls, such 
as diagonal walls, was not considered. Experimentally, 
in the Ndo.4 system Tranquada et. al. jjj report coexis- 
tence of superconducting and domain order for a range 
of dopings away from x = 0.125. For 0.05 < x < 0.12, 
the inverse domain spacing was found to vary as 2x. Be- 
yond x = 0.12 the inverse spacing remained relatively 
constant, increasing slightly as x approached 0.2. 

Here we address the question of the stability of domain 
walls at low to moderate filling, by calculating the energy 
of a domain wall in the 2D t-J model, as a function of fill- 
ing, using density matrix renormalization group (DMRG) 
H techniques. We use systems with boundary conditions 
(BCs) carefully chosen not to frustrate the domain walls. 
We also estimate the repulsive interaction between do- 
main walls, allowing us to study the domain-wall filling 
and spacing of a striped phase as a function of doping 
x. The DMRG results for energies use extrapolation to 
extract the limit of zero truncation error [Q , with up to 
1400 states per block kept. More details of the numerical 
techniques can be found in H. 



The t-J Hamiltonian in the subspacc of no doubly 
occupied sites is given by 



H = -* E ( c U-. + h - c -) + ^E( s ^ 
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Here (ij) are near-neighbor sites, s is a spin index, Si and 
cj are electron spin and creation operators, and rii — 
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c^,Cj|. The near-neighbor hopping and exchange 
interactions are t and J. We measure energies in units of 
t. We consider only J/t = 0.35 here. 

First, we consider the energetics of a single domain 
wall. Imagine a single long domain wall in the form of 
a closed loop, with a fixed number of holes. The loop 
would be a large rectangle, if domain walls prefer to be 
oriented in the (1,0) or (0,1) directions. If we assume 
that the domain wall is stable against evaporation into 
holes or pairs, then the loop will adjust its size in order 
to minimize its energy. The linear hole density pg of the 
domain wall will have an optimal value pg. We expect 
that at very low doing, any domain wall will have doping 
pg. At higher doping, repulsion between domain walls 
could lead to increased values of pg. 

Let e(pg ) be the energy per hole of a domain wall with 
density pg. Then pg minimizes e(pg). In order to measure 
e{pg), some care is needed. The important point is that a 
loop can shrink or contract without inducing frustration 
in the antiferromagnetic region inside or out. Therefore, 
in order to reduce finite size effects, we define e(pg) using 
two systems, each without frustration. The first system, 
with energy E^-w, has a domain wall and BCs favoring the 
7r-phase shift of the antiferromagnetism induced by the 
domain wall, while the second, with energy Eq, has no 
holes at all and BCs favoring no 7r-phase shift, but which 
are otherwise identical. Defining N^w a s the number of 
holes in the domain wall, we have 



Edw — E a 
A^riw 



(2) 



In practice, we use open BCs, with staggered fields on 
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the edges on either side of the domain wall to induce the 
desired antiferromagnetic order. 

-1.55 r- 



-1.60 



G— ©Domain Wall (16x6) 
1 Hole (8x8) 

1 Pair (8x8) 

■ — ■ Diagonal Domain Wall 
% — X-w/o phase shift (16x6) 
ODomain Wall (16x7) 



Q. 



-1.65 



-1.70 



-1.75 




0.0 



0.5. 



1.0 




10 20 30 40 
/ 



domain wall. Since our BCs for the doped system require 
this phase shift, for pe < 0.3 the energy per hole is quite 
high. In this case BCs without the tt phase shift give a 
lower energy, as shown by the stars. However, here we 
find that the holes bind into isolated pairs rather than a 
stripe (not shown). The two-hole energy without the 7r 
phase shift is higher than the two-hole 8x8 line because 
of larger finite-size effects on the 16 x 6 system. 

The concave nature of the 16 x 6 domain wall energy for 
0.5 < pg < 1 suggests that in this region, a long domain 
wall will phase separate into a region with p/_ = 1 and 
a region with pe = 0.5. We have directly observed this 
phase separation in a long 40 x 6 domain wall. Fig. 1(b) 
shows the density profile along the wall for pe — 0.75. 
Here the holes have separated into regions with pe = 1 
and pe ~ 0.5 — 0.6. 

In Fig. 1(a), we also show results for the energy per 
hole of a diagonal domain wall, using a tilted 12 x 7 sys- 
tem, which includes seven adjacent (1,1) lines of sites. 
Similar staggered fields for doped and undoped systems 
were applied as for the (1,0) walls shown in Fig. 1(a). 
Here the linear hole density is defined as the number 
of holes per \/2 lattice spacings, so that pe — 1 corre- 
sponds to a filled diagonal wall. The energy near pe — 1 
is slightly less than the energy of the (1,0) wall; however, 
the width of the diagonal system is slightly greater. The 
result shown for a (1,0) domain wall on a 16 x 7 system 
shows that the energies are actually nearly degenerate at 
pi = 1. At lower values of pg, (1,0) walls are lower in 
energy. 



FIG. 1. (a) Energy per hole of various hole configurations, 
as discussed in the text, (b) Linear density profile for a 40 x 6 
wall, showing phase separation. 

In Fig. 1(a) we show e(pe) as a function of pe, measured 
on a 16 x 6 system with open BCs, with the domain wall 
parallel to the x-axis, and with staggered magnetic fields 
of magnitude 0.1 applied to the top and bottom rows of 
sites. We see a minimum at pg ~ 0.5. On the same plot, 
we show by horizontal lines the energy per hole of one 
hole and of two holes placed in an open 8x8 system, 
with a staggered field of magnitude 0.1 on all four sides. 
The fact that the energy per hole of two holes is lower 
than that of a single hole indicates that two holes pair- 
bind, in agreement with exact diagonalization studies [^). 
However, the energy of a domain wall at pe is even lower, 
indicating that domain walls are stable at arbitrarily low 
doping. To judge the effects of the finite width of the 
system, results are shown for a 16 x 7 system at two 
values of pe- On these larger systems, the energy of the 
domain wall is lower [0]. 

The pe = 1 domain wall repels additional holes, leading 
to the rapid increase in e(pe) for pg > 1. For pe < 0.3, the 
holes are too far apart to induce the tt phase shift of a 
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FIG. 2. (a)-(c) Hole density, showing domain walls in a 
14x7 system. The figures are rotated by 90° for compactness, 
(d) Hole and spin densities showing a looped wall on a 20 x 8 
system with 16 holes. The diameter of the gray holes and the 
length of the arrows are proportional to 1 — (jii) and (Sf), 
respectively, according to the scales shown. 
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In Fig. 2(a)-(c) we show a system which allows either 
a (1,0) or (1,1) filled domain wall: a 14 x 7 system with 
cylindrical BCs (periodic in y, open in x), doped with 14 
holes. Here a single domain wall extends the length of the 
system. Staggered edge fields on the left and right edges 
and a small local potential on one site on each edge were 
used to pin the ends of the domain wall at specified sites. 
The initial DMRG buildup of the lattice was arranged to 
initially force an approximate domain wall of the spec- 
ified form to appear; however, later sweeps allowed the 
system to relax, although the ends of the walls remained 
pinned. Within error bars, all three domain wall config- 
urations shown in Fig. 2(a-c) have identical energies. In 
Fig. 2(a), a (1,1) domain wall wraps around the system. 
In Fig. 2(b), a wall with both (1,0) and (1,1) parts is 
present. Notice that the wall resists being situated at an 
intermediate angle. In Fig. 2(c), a (1,0) domain wall is 
present. The degeneracy of these states indicates that 
filled domain walls have the same energy, whether they 
run in the (1,0) or (1,1) directions. This suggests that 
filled domain walls might readily fluctuate or form static 
disordered configurations. 

Note that at low doping, in a strictly 2D system, one 
would expect infinite, straight (1,0) domain walls with 
filling given by the optimal filling pi w 0.5. However, in 
a system with weak coupling to other planes, in order to 
maintain long range antiferromagnetic order, it is more 
likely that domain walls would form closed loops, so that 
most of each plane would be in the dominant antiferro- 
magnetic domain. In Fig. 2(d) we show a loop on a 
20 x 8 system. In this system, staggered fields without 
any phase shifts were applied to all four sides, preventing 
a domain wall from ending on a side. Under these cir- 
cumstances, a loop forms. In a system of weakly coupled 
planes, the size of a typical loop would be set to balance 
repulsion between opposite sides of the loop and the ex- 
change cost from the coupling to adjacent planes. (If the 
exchange coupling between planes were J' ~ 10~ 5 J, and 
assuming the repulsion is given by Eq. (3), the loop size 
would be about 15-20 lattice spacings.) At higher dop- 
ings of just a few percent, the interactions between walls 
would favor the more closely packed striped phases. 

In order to understand this interaction between do- 
main walls, we have studied an L x 6 system with cylin- 
drical BCs. In this system, transverse domain walls with 
four holes wrap around the system, and are stable at low 
to moderate doping. These walls have pi = 2/3. We 
studied systems with eight or twelve holes, forming two 
or three domain walls, and studied various lengths L to 
determine the interaction between the walls. The energy 
per hole as a function of the domain wall spacing d is 
shown in Fig. 3(a). The walls repel, rather strongly at 
short distances. The solid curve in Fig. 3(a) is a simple 
exponential fit 

e(d) = e(oo) + V(d) = e(oo) + Ae'^ (3) 



with e(oo) = -1.79, A = 0.87 and w = 1.8. The source of 
the repulsion appears to be the finite width of the walls: 
the hole density distribution spreads out over several lat- 
tice spacings. In the insert to Fig. 3(a), we show the hole 
density per site as a function of l x for a 13 x 6 system with 
a single domain wall in the center. This gives the density 
profile of a wall. An isolated wall, far from boundaries in 
the L x 6 system with pi = 2/3, is site-centered. How- 
ever, walls near boundaries can be more bond-centered; 
there appears to be little energy difference between walls 
which are site centered, bond centered, or in between. 
Notice the substantial width of the wall; only 30% of the 
hole density is on the center leg. The effective mass of 
the wall as a whole seems to be very high, so that it is ef- 
fectively pinned by truncation errors in DMRG. In other 
words, we believe very little of the apparent width shown 
is due to uniform motion of the entire wall |M| . 
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FIG. 3. (a) Energy per hole of (0,1) domain walls on a 
16 x 6 system, as a function of domain wall separation. The 
inset shows the transverse hole density in one of the walls, 
(b) The energy of an array of domain walls as a function of 
doping x. 

Using the results of Fig. 1 for the energy per hole in 
a domain wall and Eq.(3) for the repulsion per hole be- 
tween domain walls, we consider the relationship between 
the domain wall spacing d and the doping x. The energy 
per site of an array of domain walls is given by 

E(p e ,d) = x[e(p e ) + V(d)}. (4) 

For pt = 1/2, d = l/(2x), while for p t = 1, d = l/x. For 
fixed pe, we define 
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E(x) = E(p e , d(x)) - x[2e(l) - e(0.5)]. 



(5) 



The latter term in Eq. (5) is like a shift in the chemical 
potential, which allows the curvature in the energy to be 
seen more easily. In Fig. 3(b), we plot E(x) for pi = 1/2 
and pi = 1 arrays of walls. Clearly, at low values of x, 
all the walls have pg = 0.5 and therefore dr 1 — 2x. For 
large values of x, all the walls have pg = 1 and therefore 



d- 1 = 



At intermediate values of x, 0.11 < x < 0.16, 



using a Maxwell construction one finds a mixture of walls 
with two different spacings di/2 and d\. However, the 
precise values of d\/2 and d\ are sensitive to V(d), which 
was determined only roughly using walls on an L x 6 
system with pg = 2/3. 
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FIG. 4. (a) Domain walls with pg = 1 and pi = 0.5 on a 
28 x 8 system, with 32 holes and cylindrical BCs. (b) Shift 
in the magnetic structure factor peak from (n, n) as a func- 
tion of x. The solid circles are experimental results[6] for 
La2-a;Sr 2; Cu04. 

To determine d\/2 and d\ more precisely, we have sim- 
ulated a 28 x 8 system, with cylindrical BCs, with four 
(0,1) pe = 0.5 walls and two (0,1) pg — 1.0 walls, and 
x = 0.14. The results are shown in Fig. 4(a). Here, the 
wall spacings naturally adjust to d\/2 and d\. We find 
di/2 = 4 and d\ = 6, implying that domain wall phase 
separation occurs for 0.125 < x < 0.17. (Note that in 
this case, the system cannot continuously adjust pg away 
from 1/2 or 1, and it may be that interactions between 
walls would shift the pi = 0.5 wall to a somewhat differ- 
ent filling.) The resulting inverse domain wall spacing is 
plotted in Fig. 4(b). Also shown are experimental results 
H for the shift in the magnetic peak in I^-^Sr^CuO^ 
Here we assume that for x < 0.03 the system has loops 
rather than stripes so that the magnetic peak would re- 
main at (tt, tt). For x > 0.125, walls with pi = 1 and 
spacing d\ give rise to an additional peak, which has not 
been observed thus far. However, as is obvious from Fig. 
4(a), the pg = 1 walls seem to be much more subject 



to disorder and fluctuations, which would significantly 
weaken and broaden this peak. Observation of this sec- 
ond peak would lend strong support to our results, as 
would a nonmonotonic shift in the peak near x > 0.2. 

In our DMRG calculations, large scale fluctuations are 
difficult to observe, so that a slowly fluctuating striped 
phase would tend to appear static to us. While this is 
certainly a disadvantage of DMRG, it makes it easier to 
rule out uniform phases — phases in which there are no 
signs of domain walls. In numerous calculations with 
various BCs, and x < 0.3, we have not seen a uniform 
phase H . We have observed apparently disordered walls, 
particularly with x = 0.15 — 0.20, but it is difficult to de- 
termine in these cases whether an ordered phase is being 
frustrated by BCs. Note that even if the pg = 1 walls are 
fluctuating or disordered, this probably does not strongly 
reduce the energy per site relative to the slightly disor- 
dered walls shown in Fig. 4(a), so that the phase separa- 
tion into regions with pg = 0.5 walls would be unaffected. 

In summary, we find that the (1,0) domain walls 
formed in the doped t-J model have a linear filling 
pg ~ 0.5 and an inverse spacing <i _1 = 2x for doping 



0.125. For x > 0.17, the domain walls have 



1 



and d -1 = x. It is tempting to identify the underdoped 
regime with the pg = 0.5 walls and the overdoped regime 
with pg = 1. We also find that pg ~ 1 (1,0) walls are 
nearly degenerate with diagonal (1,1) domain walls, sug- 
gesting that these walls may have large fluctuations. 
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